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Abstract 

The Lindblad master equation for an open quantum system with a Hamilto- 
nian containing an arbitrary potential is written as an equation for the Wigner 
distribution function in the phase space representation. The time derivative of 
this function is given by a sum of three parts: the classical one, the quantum cor- 
rections and the contribution due to the opening of the system. In the particular 
case of a harmonic oscillator, quantum corrections do not exist. 

1 Introduction 

In the last two decades, more and more interest has arisen about the problem of dissipa- 
tion in quantum mechanics and the search for a consistent description of open quantum 
systems [1-4]. Because dissipative processes imply irreversibility and, therefore, a pre- 
ferred direction in time, it is generally thought that quantum dynamical semigroups 
are the basic tools to introduce dissipation in quantum mechanics. A general form of 
the generators of such semigroups was given by Lindblad [5-7]. This formalism has 
been studied for the case of damped harmonic oscillators [8,9] and applied to various 
physical phenomena, for instance, to the damping of collective modes in deep inelastic 
collisions in nuclear physics [10-12] and to the interaction of a two-level atom with a 
single mode of the electromagnetic field [13] (for a recent review see Ref. [14]). 



The quasiprobability distribution functions have proven to be of great use in the 
study of quantum mechanical systems. They are useful not only as calculational tools, 
but can also provide insights into the connection between classical and quantum me- 
chanics. The first of these quasiprobability distributions was introduced by Wigner 
[15] to study quantum corrections to classical statistical mechanics. The Wigner dis- 
tribution function has found many applications primarily in statistical mechanics or 
in purely quantum mechanical problems, but also in areas such as quantum chemistry 
and quantum optics, collisions [16], quantum chaos [17,18] quantum fluid dynamics 
and for some aspects of density functional theory [19]. Quantum optics has given rise 
to a number of quasiprobability distributions, the most well-known being the Glauber 
P distribution [20,21]. More recently extensive use has been made of the generalized 
P distributions [22,23]. 

In the previous papers [24,25] the applicability of quasiprobability distributions to 
the Lindblad theory was explored. In order to have a simple formalism we studied the 
master equation of the one-dimensional damped harmonic oscillator as example of an 
open quantum system. This equation was transformed into Fokker-Planck type equa- 
tions for c-number quasiprobability distributions associated with the density operator 
and a comparative study of the Glauber, antinormal ordering and Wigner representa- 
tions was made. We then solved these Fokker-Planck equations, subject to different 
types of initial conditions. 

The increased activity in heavy ion physics has led to a greater interest in phase 
space distributions. Currently models as the intranuclear cascade model, "hot spot" 
model and hydrodynamical models, which emphasize the collective and transport be- 
haviour, are especially suitable for an examination in this framework [16]. We could 
also mention the nuclear reaction theory formulated in the language of the Wigner 
distribution function [26], the use of a modified Vlasov equation for the description of 
intermediate energy or high energy heavy ion collisions [27-29] , the quantum molecular 
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dynamics approach to investigate the fragment formation and the nuclear equation of 
state in heavy ion collisions [30] and the transport phenomena in dissipative heavy ion 
collisions [31]. 

The aim of the present paper is to use the quantum mechanical phase space Wigner 
distribution for a one-dimensional particle in a general potential in the framework of 
the Lindblad theory for open quantum systems. Since the Wigner distribution is a 
quantum generalization of the Boltzmann one-particle distribution, we can obtain an 
analogue of the Boltzmann- Vlasov equation for open quantum systems. 

The content of this paper is arranged as follows. In Sec. 2 we review very briefly 
the Lindblad model for the theory of open quantum systems, namely we formulate 
the basic equation of motion for the density operator. In Sec. 3 we discuss the 
quantum mechanical phase space Wigner distribution and its properties, which are 
quite useful in applications. In Sec. 4 we derive a Boltzmann- Vlasov type equation 
from the Lindblad master equation for the time dependence of the Wigner distribution. 
Finally, we discuss and summarize our results in Sec. 5. 

2 Quantum mechanical Markovian master equation 

A possibility for an irreversible behaviour in a finite system is to avoid the unitary 
time development and to consider non-Hamiltonian systems. If S is a limited set of 
macroscopic degrees of freedom and R the set of non-explicitly described degrees of 
freedom of a large system S + R, the simplest dynamics for S which describes an 
irreversible process is a semigroup of transformations introducing a preferred direction 
in time. 

In Lindblad's axiomatic formalism of introducing dissipation in quantum mechan- 
ics, the usual von Neumann-Liouville equation ruling the time evolution of closed 
quantum systems is replaced by the following Markovian master equation for the den- 
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sity operator p{t) in the Schrodinger picture [5]: 

Here, $ t denotes the dynamical semigroup describing the irreversible time evolution of 
the open system in the Schrodinger representation and L is the infinitesimal generator 
of the dynamical semigroup $ t . Using the Lindblad theorem which gives the most 
general form of a bounded, completely dissipative generator L, we obtain the explicit 
form of the most general quantum mechanical master equation of Markovian type: 

^ = ~[H,p(t)} + ^£(feW,V?] + />(*)*?])■ (2-2) 
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Here H is the Hamiltonian operator of the system and Vj, Vj are bounded operators on 
the Hilbert space 7i of the Hamiltonian. We should like to mention that the Markovian 
master equations found in the literature are of this form after some rearrangement 
of terms, even for unbounded generators. We make the basic assumption that the 
general form (2.2) of the master equation with a bounded generator is also valid for 
an unbounded generator. Since we study the one-dimensional problem, we consider 
that the operators H, Vj, Vj are functions of the basic observables p and q of the one- 
dimensional quantum mechanical system (with [q,p] = ihl, where / is the identity 
operator on H). For simplicity we consider that the operators Vj, V- are first degree 
polynomials in p and q. Since in the linear space of the first degree polynomials in p 
and q the operators p and q give a basis, there exist only two C-linear independent 
operators V\, V 2 which can be written in the form of 

Vj = ajp + bjq, J = 1,2, (2.3) 

with a,j, bj complex numbers [6]. The constant term is omitted because its contribution 
to the generator L is equivalent to terms in H linear in p and q which for simplicity 
are chosen to be zero, so that H is taken of the form 

H = H + !±(pq + qp), H = + U (q) , (2.4) 



where U(q) is the potential energy and m is the mass of the particle. With these 
choices and with the notations 

D w = | 2 M 2 > A>p = | S N*' A>? = D 9P = -| Re H a*6 i? A = -Im £ a*6 j? 

Z j=l,2 Z j=l,2 Z j=l,2 j=l,2 

(2.5) 

where a* and 6* denote the complex conjugate of dj, and b* respectively, the master 
equation (2.2) takes the following form [8,14]: 

dp i *■ i(\ — a) i(X + a) r 

* = " p| + s» ^ p " + 9P| " s» pp + pp] 

-^14, 14,?]} - frlP. \p,i>\\ + ^FIP. 14, p}}- (2.6) 

Here the quantum mechanical diffusion coefficients D pp , D qq , D pq and the friction con- 
stant A satisfy the following fundamental constraints [8,14]: 

i) D pp > 0, ii) D qq > 0, iii) D pp D qq - D pq 2 > X 2 h 2 /A. (2.7) 

The necessary and sufficient condition for L to be translationally invariant is p — A 
[6,8,14]. In the following general values for A and p will be considered. 

3 Wigner distribution function 

In the phase space picture, the quantum corrections become transparent and a smooth 
transition from quantum to classical physics is encountered. This picture is particularly 
suitable for obtaining quantum mechanical results in situations where a good initial 
approximation comes from the classical result and also for deriving classical limits of 
quantal processes. 

A quantum mechanical particle is described by a density matrix p and the average 
of a function of the position and momentum operators A(p, q) (we consider the one- 
dimensional case) is 

< A >= Tr(ip), (3.1) 



or, in terms of a quasiprobability distribution 4>(p,q), 

<A>= J dp J dqA{p, q)(f>(p, q), (3.2) 

where the function A(p, q) can be derived from the operator A(p, q) by a well-defined 
correspondence rule [32]. This allows one to cast quantum mechanical results into a 
form in which they resemble classical ones. In the case where <f> in Eq. (3.2) is chosen 
to be the Wigner distribution, then the correspondence between A(p, q) and A is that 
proposed by Weyl [33], as was first demonstrated by Moyal [34]. 

The Wigner distribution function, which represents the Weyl transform of the 
density operator, is defined through the partial Fourier transform of the off-diagonal 
elements of the density matrix: 

W(p, q) = ^-Jdy<q- y\p\q + y> e 2 ™' h (3.3) 

and satisfies the following properties: 

i) W(p, q) is real, but cannot be everywhere positive, 

n)fdpW(p,q)=<q\p\q>, (3.4) 
fdqW(p,q) =<p\ P \p>, (3.5) 
fdpfdqW(p,q)=Trp=l, (3.6) 

iii) W(p, q) is Galilei invariant and invariant with respect to space and time reflec- 
tions, 

iv) In the force free case the equation of motion is the classical (Liouville) one 

dW _pdW 

dt ~ mdq- [6J) 
Wigner has shown that any real distribution function as long as it satisfies properties 
ii) assumes also negative values for some p and q. 

The classical phase space function A(p, q) corresponding to the quantum operator 
A is given by the so-called Wigner-Moyal transform: 

A(p, q)= [ dze tpz/h <q- \z\A\q + \z >, (3.8) 



so that / / dpdqA(p, q) = 2tc%TiA. Clearly, Eq. (3.3) is a special case of Eq. (3.8) 
for the density operator, i.e. A = p and W(p, q) is the phase space function which 
corresponds to the operator p/2nh. The function A(p,q) is known as the Wigner 
equivalent of the operator A. The relation which expresses the Wigner equivalent F of 
the product of operators F = AB in terms of the Wigner equivalents of the individual 
operators A and B is the following: 

F(p,q) =A(p,q)(exp^)(p,q) = B(p,q)(exp(~))A(p,q), (3.9) 

where A (essentially the Poisson bracket operator) is given by 

dp dq dq dp 

and the arrows indicate in which direction the derivatives act. From Eq. (3.9) the 
Wigner equivalent of the commutator [A, B] follows directly: 

([A,B])(p,q) = -2iA(p,q)(sm^)B(p,q). (3.11) 

Next we recall the time dependence of the Wigner distribution function W(p, q,t) for 
a closed system [32]. Instead of the von Neumann-Liouville master equation we have 
the following quantum Liouville equation which determines the time evolution of the 
Wigner phase space distribution function (Wigner equivalent of the density operator): 

tft^ = H{p, g)(exp ^)W(p, q) - W(p, g)(exp ^)H(p, q), (3.12) 



or, cf. Eq. (3.11), 



= -2H(p,q)(sm^)W(p,q), (3.13) 



where H(p, q) is the Wigner equivalent of the Hamiltonian operator H of the system. 
Note that if we take the h — > limit of this equation, we obtain the classical Liouville 
equation (Vlasov equation). 
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4 The equation of motion for the Wigner distribu- 
tion of open systems 

In two previous papers [8,24] the evolution equations for the Wigner function were 
obtained for the damped harmonic oscillator within the Lindblad theory for open 
quantum systems. In this Section we obtain the equation for the Wigner function of 
a one-dimensional system with an arbitrary potential. In other words, we obtain, by 
means of the Wigner function, the phase space representation of the general Lindblad 
master equation. 

The time evolution of the Wigner distribution function corresponding to the Lind- 
blad master equation (2.2), can be obtained from Eq. (3.12) by adding in the right- 
hand side the Wigner equivalent corresponding to the opening part of Eq. (2.2), i.e. 
the sum of commutators. By using formulas (3.9) and (3.11), we obtain the following 
evolution equation for the Wigner distribution: 

dW 2 ft A 
= --H(sm—)W 



dt h y 2 



1 \- 

3 



^)W(exp ^ )V * - V7(exp ^(exp - W(e W ^/(exp 

(4.1) 



where Vj and V* are the Wigner equivalents of the operators Vj and Vj , respectively, 
and V* is the complex conjugate of Vj. If the operators Vj are taken of the form (2.3), 
then Eq. (4.1) becomes the evolution equation for the Wigner function corresponding 
to the master equation (2.6): 

dW 2 HA d d 

— = --„ {sinT) w + x rq ( <1 w) + x^( P w) 

„ d 2 W „ d 2 W ^ d 2 W , AnS 
+D„ W + D„— + 2D M — . (4.2) 

Here it is easily to remark that the first term on the right-hand side generates the 
evolution in phase space of a closed system and gives the Poisson bracket and the higher 
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derivatives containing the quantum contribution, and the following terms represent the 
contribution from the opening (interaction with the environment). 

For a Hamiltonian operator of the form (2.4) and by assuming a Taylor expansion 
of the potential U, i.e. U (q) is an analytic function, this equation takes the form: 

dW p dW OUdW ~ lXn (h) 2n d 2n+1 U(q) d 2n+1 W „ TT , AnS 

= — — 1 h > (— 1) — — h CW A3) 

dt m 8q dq dp ^ ' 2 2 «(2n+l)! dq 2n+1 dp 2n+l ' 1 ' ' 

where we have introduced the notation 

d d d 2 W d 2 W d 2 W 

C = (A - ^)-(qW) + (A + »)-( P W) + D qq — + D pp — + 2D pq —. (4.4) 

If Eq. (4.3) had only the first two terms on the right-hand side, W would evolve along 
the classical flow in phase-space. The terms containing A and \i are the dissipative 
terms. They modify the flow of W and cause a contraction of each volume element 
in phase space. The terms containing D pp , D qq and D pq are the diffusive terms and 
produce an expansion of the volume elements. The diffusion terms are responsible 
for noise and also for the destruction of interference, by erasing the structure of the 
Wigner function on small scales. The sum term (the power series), together with the 
first two terms make up the unitary part of the evolution. Hence, up to corrections 
of order h 2 , unitary evolution corresponds to approximately classical evolution of the 
Wigner function. It is partly for this reason that the evolution of a quantum system 
is most conveniently undertaken in the Wigner representation. The higher corrections 
can often be assumed as negligible and give structures on small scales. There are, 
however, important examples where they cannot be neglected, e.g., in chaotic systems 
(see the discussion below). From Eq. (4.3) it is clear that, as a consequence of the 
quantum correction terms with higher derivatives, the Wigner function of a non-linear 
system does not follow the classical Liouville flow. The higher derivative terms are 
generated by the nonlinearities in the potential U(q). 

The question of the classical limit of the Wigner equation is of more than formal 

9 



interest [35]. There are two well-known limits in which Eq. (4.3) can go over into a 
classical equation: 1) when U is at most quadratic in q and 2) when h — > 0. 

Because of the extra diffusion terms, we get yet a third classical limit: In the 
limit of large D pp , the diffusive smoothing becomes so effective that it damps out all 
the momentum-derivatives in the infinite sum and Eq. (4.3) approaches the Liouville 
equation with diffusion, an equation of Fokker-Planck type. This is an example of how 
macroscopic objects start to behave classically (decoherence), since the diffusion coeffi- 
cients are roughly proportional to the size of these objects. Thus an object will evolve 
according to classical dynamics if it has a strong interaction with its environment. 
The connection between decoherence and transition from quantum to classical in the 
framework of the Lindblad theory for open quantum systems will form the subject of 
another paper [36]. 

In the following we shall consider Eq. (4.3) for some particular cases. 

1) In the case of a free particle, i.e. U(q) = 0, Eq. (4.3) takes the form: 

dW v dW 
at m oq 

2) In the case of a linear potential U = ^q (where for example 7 = mg for the free 
fall or 7 = eE for the motion in a uniform electric field), one gets 

dW p dW dW , , 

— = -^— + 1 — + CW. 4.6 
at m oq op 

3) In the case of the harmonic oscillator with U = mu 2 q 2 /2, Eq. (4.3) takes the 
form: 

dW pdW 2 dW _ , n . 

— = + mu 2 q—- + CW. 4.7 

at m oq Op 

An analogous equation can be obtained in the case of the motion on an inverted 
parabolic potential U(q) = —mK 2 q 2 /2. Since the drift coefficients are linear in the 
variables p and q and the diffusion coefficients are constant with respect to p and 
q, Eqs. (4.5)-(4.7) describe an Ornstein-Uhlenbeck process [37-39]. Eqs. (4.5)-(4.7) 
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are exactly equations of the Fokker-Planck type. Eq. (4.7) was extensively studied 
is the literature [8,14,25] and it represents an exactly solvable model. It should be 
stressed that not every function W(p, q, 0) on the phase-space is the Wigner transform 
of a density operator. Hence, the quantum mechanics appears now in the restrictions 
imposed on the initial condition W(p, q, 0) for Eq. (4.3). The most frequently used 
choice for W(p, q, 0) is a Gaussian function and Eqs. (4. 5)- (4. 7) preserve this Gaussian 
type, i.e., W(p,q,t) is always a Gaussian function in time, so that the differences 
between quantum and classical mechanics are completely lost in this representation of 
the master equation. 

4) In the case of an exponential potential U(q) = aexp(—j3q), the Wigner equation 
is an infinite order partial differential equation 



(4.8) 

but in the case of a potential of the finite polynomial form U (q) = Y,n=i a nq n , the sum 
keeps only a certain number of derivate terms. As an illustration of this remark, we 
consider an anharmonic oscillator with the potential U an h(q) = mu 2 q 2 /2 + Cq A . In this 
case the Wigner equation (4.3) becomes 

dW pdW . 2 3 ,dW d 3 W „ Jir 

— = + ijnu 2 q + 4C7g 3 )-— - Ch 2 q—- + CW. 4.9 

at m oq op op 6 

The above equation has one term with third derivative, associated with the nonlinear 
potential U anh . In fact, the first three terms on the right-hand side of Eq. (4.9) give 
the usual Wigner equation of an isolated anharmonic oscillator. The third derivative 
term is of order h 2 and is the quantum correction. In the classical limit, when this 
term is neglected, the Wigner equation becomes one of the Fokker-Planck type. 
5) For a periodic potential U(q) = U cos(kq) : 
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d 2n+l U 



(-l) n k 



dU 
dq 



(4.10) 



Q q 2n+1 



and we obtain 



dW 



p dW dU 5W 
JL 1 

m dq dq Sp 



+ cw, 



(4.11) 



dt 



where 



SW _ W(q, p + hk/2, t) - W(q, p - hk/2, t) 



(4.12) 




Eq. (4.11) takes a simpler form when hk is large compared to the momentum spread 
Ap of the particle being considered, i.e. when the spatial extension of the wave packet 
representing the particle is large compared to the spatial period of the potential. Im- 
posing the condition hk 3> Ap on Eq. (4.12), one sees that SW/Sp is small for any p 
that yields an appreciable value of the Wigner distribution function, i.e. 5W/5p ~ 
for all practical purposes [40]. Eq. (4.11) is then reduced to the equation (4.5) for a 
free particle moving in an environment. 

From the examples l)-3) we see that for Hamiltonians at most quadratic in q and 
p, the equation of motion of the Wigner function contains only the classical part and 
the contributions from the opening of the system and obeys classical Fokker-Planck 
equations of motion (4.5)-(4.7). In general, of course, the potential U has terms of 
order higher than q 2 and one has to deal with a partial differential equation of order 
higher than two or generally of infinite order. When the potential deviates only slightly 
from the harmonic potential, one can still take the classical limit h — > in Eq. (4.3) 
as the lowest-order approximation to the quantum motion and construct higher-order 
approximations that contain quantum corrections to the classical trajectories using 
the standard perturbation technique [40]. 

Eq. (4.3) can also be rewritten in the form 



dW 



p dW dUeff dW 
m dq dq dp 



+ CW, 



(4.13) 



dt 
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where the effective potential is defined as [40] 



dU e ff dW 
dq dp 



dU dW 



+ E(-i) 



oo 



n 



(h) 2n d 2n+1 u( q ) d 2n+1 w 



(4.14) 



dq dp 



n=l 



2 2n (2n+ 1)! dq 2n+1 dp 2n+1 ' 



Then the phase-space points move under the influence of the effective potential U e ff. 
We note that Eq. (4.14) indicates that only when at least an approximate solution for 
W is known, the effective potential can be determined. If U does not deviate much 
from the harmonic potential, the zeroth-order approximation for W can be taken as 
that resulting from classical propagation. The main limitation of the effective potential 
method is that it can be applied only to systems whose behaviour is not much different 
from that of the harmonic oscillator or the free wave packet, for example a low-energy 
Morse oscillator and an almost free wave packet slightly perturbed by a potential step 
or barrier [40]. The effective potential method may possibly be applied to a collision 
system [40]. The quantum dynamics of atomic and molecular collision processes has 
been described rather successfully in the past by means of classical trajectories [16,40- 
42], i.e. the collision system is one for which an approximate solution of W can 
be considered as known and thus the perturbative scheme of the effective potential 
method can be employed. In Refs. [43-45] similar equations were written for the 
Wigner function for a class of quantum Brownian motion models consisting of a particle 
moving in a general potential and coupled to an environment of harmonic oscillators. 

Phase space provides a natural framework to study the consequences of the chaotic 
dynamics and its interplay with decoherence [46-50]. Eq. (4.3) could be applied in 
order to investigate general implications of the process of decoherence for quantum 
chaos. Since decoherence induces a transition from quantum to classical mechanics, it 
can be used to find the connection between the classical and quantum chaotic systems. 
In this case U (q) is the potential of a classically chaotic system, coupled to the external 
environment. For this purpose, a particular case of Eq. (4.3) was utilized by Zurek 
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and Paz [51] in the special case of the high temperature limit of the environment: 
/i = A = 7, where 7 is the relaxation rate, D qq = 0, D pq = and the diffusion coefficient 
D pp = D — 2m^fkBT (T is the temperature of the environment). In this model the 
symmetry between q and p is broken, and coupling with the environment through 
position gives momentum diffusion only [52]. Zurek also used a similar equation in 
which the diffusion term is symmetric, namely D(dp p + dg)W [51]. 

It is sometimes argued that the power series involving third and higher derivative 
terms may be neglected. Paz and Zurek [51] have argued that the diffusive terms 
may smooth out the Wigner function, suppressing contributions from the higher-order 
terms. When these terms can be neglected, the Wigner function evolution equation 
(4.3) then becomes 

dw pdw dudw rijr „ . 

at m oq oq op 
In the previous considered particular case of a thermal bath and if A = /i — 7, D qq = 

D pq = 0, Dpp = D — 21717 ksT, this equation becomes a Kramers equation [45,53]: 

dW _ _p_dW_ + dU dW + 2 ^d_ + 2 T ^W_ 
dt in dq dq dp dp dp 2 

It posseses the stationary solution [45] 

P 2 U(q) 



W(p, q) = N exp 



(4.17) 



2mk B T k B T 

where N is a normalization factor. As was discussed by Anastopoulos and Halliwell 
[45], this will be an admissable solution, i.e., it is a Wigner function only if the potential 
is such that exp[—U(q)/D] is normalizable. This requires U(q) — > 00 as q — > ±00 faster 
than In \ q\ . In this case the stationary distribution is the Wigner transform of a thermal 
state p = Z~ x exp(—H /kT) with Z = Tr exp(— H /kT) for large temperatures. The 
general results of [53] then show that all solutions of Eq. (4.16) approach the stationary 
solutions (4.17) as t — > 00. Hence, to the extent that Eq. (4.16) is valid, all initial 
states tend towards the thermal state in the long-time limit [45] . 
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5 Summary and outlook 

The Lindblad theory provides a selfconsistent treatment of damping as a possible ex- 
tension of quantum mechanics to open systems. In the present paper we have shown 
how the Wigner distribution can be used to solve the problem of dissipation for some 
simple systems. From the master equation we have derived the corresponding Vlasov 
or Fokker-Planck equations in the Wigner W representation. The Fokker-Planck equa- 
tions which we obtained describe an Ornstein-Uhlenbeck process. The Wigner repre- 
sentation discussed provides a classical conceptual framework for studying quantum 
phenomena and enables one to employ various methods of approximation or expansion 
used in classical cases to problems of the quantum domain. By means of the Wigner 
function one can calculate mean values of physical quantities that may be of inter- 
est, using a classical-like phase space integration, where position and momentum are 
treated as ordinary variables rather than as operators. 

The phase-space formulation of quantum mechanics represents an alternative to the 
standard wave mechanics formulation of the Lindblad equation. The main difficulty 
with the phase-space formulation is that the time development of the phase-space 
distribution (Wigner) function is given in terms of an infinite-order partial differential 
equation (see Eq. (4.3)). What matters, however, is not whether the equation can 
be solved exactly, but whether a reasonable set of approximations can be introduced 
that make it possible to obtain an approximate but still accurate solution to the 
equation. The Wigner phase-space formulation is particularly useful when the classical 
description of the system under consideration is adequate and quantum corrections to 
the classical description are only desired for higher accuracy. 

Acknowledgment. One of us (A.I.) whould like to express his sincere gratitude 
to Professor W. Scheid for the hospitality at the Institut fur Theoretische Physik in 
Giessen. 
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